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Abstract 

Using monomer-resolved Molecular Dynamics simulations and theoretical arguments based on 
the radial dependence of the osmotic pressure in the interior of a star, we systematically investigate 
the effective interactions between hard, colloidal particles and star polymers in a good solvent. The 
relevant parameters are the size ratio q between the stars and the colloids, as well as the number 
of polymeric arms / (functionality) attached to the common center of the star. By covering a wide 
range of q's ranging from zero (star against a flat wall) up to about 0.75, we establish analytical 
forms for the star-colloid interaction which are in excellent agreement with simulation results. A 
modified expression for the star-star interaction for low functionalities, / < 10 is also introduced. 
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I. INTRODUCTION 



The study of the structural and thermodynamic properties of the polymeric state of mat- 
ter has a long history in physics, which started with the pioneering work of Flory [I], [|, |3] . 
In the traditional, "polymeric approaches" to the matter, the chain nature of the macro- 
molecules involved is in the foreground. However, in the last few years, alternative, com- 
plementary considerations have emerged that can loosely be called "colloidal approaches". 
Here, one envisions the polymer chains as diffuse, spherical objects and the chain nature of 
the molecules does not explicitly appear in the formalism. Instead, in a first step, almost all 
of the monomeric degrees of freedom are thermodynamically traced out of the problem [|J. 
Thereby, the polymers are replaced either by their centers of mass or by one of the monomers 
along their backbone, typically the end- or the central monomer. In this way, effective inter- 
actions between the polymers naturally arise @], which implicitly include the effects of the 
traced-out monomers and typically have the range of the chain radius of gyration R g . The 
earliest such approach dates back to the work of Asakura and Oosawa |5], ||] , and Vrij , who 
modeled polymer chains as penetrable spheres. These models pertain mostly to Gaussian, 
i.e., ideal chains and are semi- quantitative. More systematic approaches have appeared in 
the recent years, in which self-avoiding chains are modeled and effective interactions among 
them are derived by means of simulations || or theory || . The gain from adopting such an 
alternative view is twofold: on the one hand, one has the possiblity of looking at the same 
problem from a different angle; on the other hand, tracing out the monomers reduces the 
complexity of the problem by a factor N, the degree of polymerization of the chains ||TU| . 

A physical system where the colloidal approach finds an intuitive and natural application 
is that of star polymers [jlTfl . These macromolecular entities are synthesized by covalently 
attaching / polymeric chains on a common center. In this way, hybrid particles between 
polymers and colloids can be constructed, which naturally bridge the gap between these two 
common states of soft matter. The number of arms /, also known as functionality of the stars, 
allows us to go from free chains (/ = 1,2) to stiff, spherical particles (/ ^> 1). Effective 
interactions between star polymers in good [F2| and B-solvents [|13| have been recently 
derived and the validity of the former has been confirmed through extensive comparisons 
with experiments |I2], [14], [15| and simulations |T6| , |IT| . Extensions to polydisperse stars [T3[ 
as well as to many-body forces in dense star polymer solutions [pl| have also been recently 
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carried out. 

In this work, we wish to carry these considerations one step further by looking at a 
two- component system of star polymers in good solvent conditions and hard, spherical, 
colloidal particles. Though the star-star interaction is readily available and the colloids can 
be modeled as hard spheres, the effective cross interaction between star polymers and colloids 
is still missing. It is the purpose of this work to present theoretical and simulation results and 
to furnish analytic expressions for the force and/or the effective interaction acting between 
a star polymer and a spherical, colloidal particle for a large range of size ratios between 



the two. The theoretical approach is inspired by the eariler considerations of Pincus f20 



regarding the force acting between a star and a flat wall but are made more precise here 
and they are also extended to include the effects of curvature. The rest of this paper is 
organized as follows: in section II we present the general theoretical approach, both for flat 
and curved surfaces and derive analytic expressions for the star-colloid force which include 



a handful of undetermined parameters. In section |T| we compare those with the results 
of monomer-resolved Molecular Dynamics simulations and determine the free parameters in 
order to achieve agreement between theory and simulation results. In section [IV] we present 
a modified version of the star-star potential which is valid for very low arm numbers, / < 10, 
and in section [V] we summarize and conclude. 



II. THEORY 



Let us first define the system under consideration and its relevant parameters. We con- 
sider a collection of star polymers with functionality / and hard, spherical colloidal particles, 
the interaction between the latter species being modeled through the hard sphere (HS) po- 
tential. By considering two isolated members of each species, i.e., one star and one colloid, 
our goal is to derive the effective interaction between the two. The colloids have a radius 
R c , which is a well-defined length scale. 

The stars, on the other hand, are soft, hairy balls without a sharply defined boundary 
and this leads to some freedom in defining length scales characterizing their spatial extent. 
The experimentally measurable length scale that naturally arises from small-angle neutron- 
or X-ray-scattering experiments (SANS or SAXS) is the radius of gyration R g of the stars 
and the associated diameter of gyration cr g = 2R g . For the theoretical investigations on 
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the subject, however, another length scale turns out to be more convenient, namely the so- 
called corona radius R s of the star or the associated corona diameter o~ s = 2i? s . The corona 
radius arises naturally in the blob model for the conformation of isolated stars, introduced 
by Daoud and Cotton pi| . According to the Daoud-Cotton picture, the bulk of the interior 
of a star in good solvent conditions (and for sufficiently long arm chains), consists of a region 
in which the monomer density profile c(s) follows a power-law as a function of the distance 
s from the star center, namely: 



3f£r 4/ Vv3 



c(s) ~ a | - ) v 
a 



f /3 , (1) 



with the monomer length a, the excluded volume parameter v and the reduced excluded 
volume parameter v = v/a?. Outside this scaling region, there exists a diffuse layer of 
almost freely fluctuating rest chains, in which the scaling behavior of the monomer profile 
is not any more valid. We define the corona radius R s of the star as the distance from the 
center up to which the scaling behavior of the monomer density given by Eq. (0) above holds 
true. In what follows, we define the size ratio q between the stars and the colloids as: 

^7T (2) 

tic 

In addition, the interior of the star forms a semidilute polymer solution in which scaling 
theory |22| predicts that the osmotic pressure II scales with the concentration c as 11(c) ~ 



c 9 / 4 . Combining the latter with Eq. (|T|) above, we obtain for the radial dependence of the 
osmotic pressure of the star within the scaling regime the relation: 

U(s)^k B Tf/ 2 s- 3 (s<R s ). (3) 

No relation for the osmotic pressure IT(s) for the diffuse region s > R s is known to date. 
It is indeed one of the central points of this work to introduce an accurate ansatz for the 
latter, one that will allow us also to derive closed formulas for the effective force between a 
star and a hard object. This is the subject we examine below. 



A. A star polymer and a flat wall 



We begin by examining the simplest case, in which a star center is brought within a 
distance z from a hard, flat wall, as depicted in Fig. [l|. Going back an idea put forward 
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some ten years ago by Pincus p0| , we can calculate the force F sw (z) acting between the 
polymer and the wall by integrating the normal component of the osmotic pressure II(s) 
along the area of contact between the star and the wall. In the geometry shown in Fig. [I], 
this takes the form: 

F sw (z) = 2n / n(s)costfydy. (4) 

Jy=0 

Using z = s cos$ and y = z tan^ we can transform Eq. (Q) into: 

F sw (z) = 2nz / n(s)ds. (5) 



Eq. (|5|) above implies immediately that, if the functional form for the force F sw (z) were to 
be known, then the corresponding functional form for the osmotic pressure IT(z) could be 
obtained through: 

H(z) oc -± ( . (6) 



dz \ z 

To this end, we now refer to known, exact results regarding the force acting between a 



flat wall and a single, ideal chain whose one end is held at a distance z from a flat wall p3 
There, it has been established that the force F&P(z) is given by the relation: 



FM(z) = k B T^hL[ert(l)], (7) 

where erf(x) = Ij^/H f*e~ ta dt denotes the error function and L is some length scale of 
the order of the radius of gyration of the polymer. Carrying out the derivative and setting 
erf(x) = 1 for x ^> 1, we obtain a Gaussian form for the chain- wall force at large separations: 

i^) = ^exp(--^) (z»L). (8) 

We now imagine a star composed of ideal chains. As the latter do not interact with each 
other ("ghost chains") the result of Eq. @ holds for the star as well. Going now to self- 
avoiding chains, we assert that, as the main effect giving rise to the star- wall force is the 
volume which the wall excludes to the chains, rather than the excluded volume interactions 
between the chains themselves, a relation of the form (||) must also hold for the force F sw (z) 
between a wall and a real star, but with the length scale L replaced by the radius of gyration 
or the corona radius of the latter and with an additional, /-dependent prefactor for taking 
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into account the stretching effects of the / grafted polymeric chains. From Eqs. (§) and 
it now follows that 



n( s )oc^f 3 + 4) ex P (-7?) (0) 



L 

L V^ 2 ' L 2 J^\ L 2 
The full expression for IT(s) now follows by combining Eq. (^), valid for s < R s , with Eq. 
(P), valid for s ^> L = i? s , and matching them at s = R s . The local osmotic pressure n(s) 
is the interior of a star polymer, as a function of the distance s from its center has hence 
the functional form: 

s~ 3 for s < R s ; 

U(s) = Af/ 2 k B T ~ (10) 

[ (jt + 2k 2 ) £ exp [-k 2 {s 2 - R 2 )} for a > R s , 

where A and k = L^ 1 are free parameters; it is to be expected that k = 0(R~ x ), as we will 
verify shortly On the other hand, £ must be chosen to guarantee that n(s) is continuous at 
s = R s , resulting into the value: 



Eq. ([10]) above concerns the radial distribution of the osmotic pressure of an isolated star. 
The question therefore arises, whether this functional form for the osmotic pressure can be 
used in order to calculate the force between a star and a flat wall also in situations where 
the star-wall separation is smaller than the radius of gyration of the star, in which case 
it is intuitively expected that the presence of the wall will seriously disturb the monomer 
distribution around the center and hence also the osmotic pressure. In fact, it is to be 
expected the osmotic pressure is a function of both the star-wall separation z and the radial 
distance s, whereas in what follows we are going to be using Eq. ([|) together with Eq. (pX|, 
in which Tl(s) has no ^-dependence itself. However, it turns out that this is an excellent 
approximation. On the one hand, it is physically plausible for large star-wall separations, 
where the presence of the wall has little effect on the segment density profile around the 
star center and the ensuing osmotic pressure profile. On the other hand, also at very small 
star-wall separations, the scaling form Tl(s) ~ s~ 3 continues to be valid. To corroborate this 
claim, we proceed with some arguments to this effect. 

First, we refer once more to known, exact results concerning the radial distribution of the 



pressure on a hard wall arising from an ideal chain grafted on it [24], a situation similar to 



holding one end of a chain at a distance very close to the wall surface. The pressure 111(1(3) 



reads as |24 



1 1 / s 2 + a 2 



s 2 + a 2 

Am 



(12) 



with the segment length a, indicating that in the regime a < s < i? g indeed the scaling 
Ilid(s) ~ s~ 3 holds. 

Second, we can employ a scaling argument, asserting that, on dimensional grounds, the 
osmotic pressure exerted by a star on a nearby flat wall and held at a distance z from it, must 
be of the form H(s, z) = k B TR~ 3 h(s/R g , z/R g ), with some scaling function h(x, y); univer- 
sality arguments dictate that the segment length a should not appear in the dimensional 
analysis and hence s, z and R g are the only relevant length scales for this problem. Now, 
for small star- wall separations, z <C fi g , we replace the second argument of this function 
by zero. Moreover, we assert that, as the dominant contribution to the osmotic pressure 
for distances s < R g comes from the first few monomers along the chains colliding with 
the wall, the degree of polymerization N of the chains should be irrelevant if the chains 
are long. Hence, all independence of the pressure should drop out, with the implication 
h(x, 0) ~ x~ 3 for i< 1 and hence IT(s) ~ s -3 in this regime. 

Third, we point out that bringing a star with / arms at a small distance to a flat wall, 
creates a conformation which is very similar to one of an isolated star with 2 /-arms, as 
shown in Fig. 0. Hence, it is not surprising that at small star-wall separations, one recovers 
for the radial dependence osmotic pressure the scaling laws pertinent to an isolated star. 

Finally, by inserting Eq. (|T0|) into Eq. (|5|) and carrying out the integration, we find that 
for small star-wall distances, z <C R s , the force scales as F sw (z) ~ (feTj/z, thus giving rise 
to a logarithmic effective star-wall potential V sw (z) ~ —kBTln(z/R s ). The latter is indeed in 
full agreement with predictions from scaling arguments arising in polymer theory pO, 25. 26 



This is a universal result, in the sense that it also holds for single chains, be it real or ideal, as 
it can also be read off from the exact result, Eq. (0), using the property erf (x) ~ x for x — > 0. 
Thus, the proposed functional form for the osmotic pressure, Eq. (|T0D , combined with Eq. 
(||) for the calculation of the effective force, has the following remarkable property: it yields 
the correct result both at small and at large star-wall distances and therefore appears to be 
a reliable analytical tool for the calculation of the effective force at all star-wall distances. 
At the same time, it contains two free parameters, A and k which allow some fine tuning 



when the predictions of the theory are to be compared with simulation results, as we will do 
below. Yet, we emphasize that this freedom is not unlimited: on physical grounds, k must 
be of the order of R" 1 and A must be a number of order unity for all functionalities /, as 
the dominant, / 3//2 -dependence of the osmotic pressure prefactor has been already explicitly 
taken into account in Eq. (]T0|) . 

We are now in a position to write down the full expression for the star-wall force, by 
using Eqs. (||) and ([U]). The result reads as: 

RsFsUz) _ A/3 / 2 h + j(^-l) ior z<R s ; ^ 

kBT ^exp[-K 2 (z 2 -R 2 s )} for z > R s . 



Note the dominant, ~ l/z-dependence for z — > 0. Accordingly, the effective interaction 
potential V sw (z) between a star and a flat, hard wall held at a center-to-surface distance z 
from each other reads as: 

^U^A/-j^ n(i) - ( *- 1)K " )+C (M) 
I ( eric(K,z) / eric(K,R s ) for z > R s , 

with the inverse temperature (3 = (ksT) -1 , the additional constant 

( = exp( K 2 J R s 2 ) erfcMs) (15) 
kR s 

and the complementary error function erfc(x) = 1 — erf(x). This completes our theoretical 
analysis of the star polymer-wall force and the ensuing effective interaction potential. The 
comparison with simulation data and the determination of the free parameters in the theory 
will be discussed in section [TTT| . We now proceed with the calculation of the effective force 
between a star and a spherical hard particle, where effects of the colloid curvature become 
important. 



B. A star polymer and a spherical colloid 

We apply the same idea as for the case of the hard wall: the effective force acting at the 
center of the objects is obtained by integrating the osmotic pressure exerted by the polymer 
on the surface of the colloid. In Fig. [3|, the geometrical situation is displayed: within the 
corona radius of the star polymer R s = cr s /2, the osmotic pressure is determined by scaling 
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laws; the outer regime is shadowed and signifies the Gaussian decay of the osmotic pressure. 
At center-to-surface distance z (center-to-center distance r = z + R c ), the integration of the 
osmotic pressure is carried out over the contact surface between star and colloid. Taking 
into account the symmetry of the problem, e.g., its independence of the azimuthal angle, we 
obtain the force F sc (z) between the star and the colloid as: 



'max 



F sc {z) = 2tcR 2 c / d9 sin#IT(s)costf, (16) 
Jo 

where $ and 9 are polar angles emanating from the center of the star polymer and the 
colloid, respectively. The variables d and 9 can be eliminated in favor of the variable s, which 
denotes the distance between the center of the star and an arbitrary point on the surface 
of the colloid. This elimination is achieved by taking into consideration the geometrical 
relations (see Fig. [3]): 

s sin = R c sin 9 (17) 

and 

s cos$ + R c cos6 = R c + z. (18) 

Eqs. (|16|) ) (|17D and ( |18D yield for the star-colloid effective force the transformed integral: 

F 8c (z) = J [(z + R c f - Rl + s 2 } U(s) (19) 

The maximum integration distance, s max , depends geometrically on 9 max , as well as on the 
distance z of the star polymer to the surface of the colloid and on R c . The relation reads as 



[z + R C (1- cos 6» max )f + (Rc sin 6> max ) 2 



= -y/[qz + R g (l-cose max )Y+ (R g sm9 max ) 2 . (20) 

By introducing Eq. ( ]TD| ) into Eq. (|T^p, an analytic expression for the effective force follows, 
which reads as 



\2 



F sc (z) Af/ 2 R C J [{z + R c f - R 2 C ] - ^ + *i(i? s )J - ln(f ) + ^(Rs) for z < R s] 

[(z + R c ) 2 - R 2 C ] Vi(z) + V 2 (z) for z > R s . 

(21) 



k B T ~ (z + R ( 
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Here, the functions *&i(x) and ^O^) are given by: 



*!(x) = |-exp {k 2 R 2 s ) 



1 / 2 2\ 1 / 2 2 \ 

-exp (-« x ) exp {-k s max ) 

A "J max 



(22) 



and 



*2(a0 = ^exp (k 2 R 2 s ) 



[erf (Ks max ) - erf (kx)] + xexp (-k 2 x 2 ) - s max exp (-/« 2 s max ) 



/v 



(23) 



where £ is given by Eq. fllTI). Note that, for small distances, both regimes of the osmotic 
pressure contribute to the integral, whereas for larger distances, z > R s , only the Gaussian 
decay does so. Due to the additional dependence of s max on the distance z, [see Eq. (0)], an 
analytical expression for the effective potential V BC (z), analogous to Eq. (|T4l ) for the flat-wall 
case, is not possible here. 

Some remarks regarding F sc (z) are necessary. First, for small separations z, the force 
scales as F sc (z) ~ {ksT)/z^ the same behavior found for the flat-wall case. Once more, we 
obtain the universal result mentioned above, which has been shown to be also valid for an 
ideal chain whose one end is held at a distance z from the surface of a hard sphere. Indeed, 
for this case the force is given by the exact relation [ 27| : 



k B T-^- In 
oz 



1 - 



Rr 



z + R c 



erfc 



L 

(id), 



(24) 



(k B T)/zioiz -> 0. 



with L being a length scale of order R g . Eq. (|24"D above, yields F s 

Second, let us consider the limit of small size ratios q = R g /R c . As can be seen from 
Eq. (^), the upper integration limit s max scales as R g /q, whereas the decay parameter n 
is of the order Rg 1 - Hence, Ks max ~ q~ x , with the implication that for small enough g's, 
the argument ns max in the error function and in the Gaussian in Eqs. ( p2|) and (|23|) can be 
replaced by infinity. As we will shortly see, this is an excellent approximation up to q < 0.3, 
as both erf(x) and exp(— x 2 ) approach their asymptotic values for x —>■ oo rapidly. Then, 
the implicit ^-dependence of the force F BC (z) through s max drops out and a ^-integration of 
the latter can be analytically carried out to obtain an effective star polymer-colloid potential 
V°°(z) which reads as 



A/ 



3/2 



Rr 



z + R c 



-ln(i)-(^-l)(e-|)+C hrz<R s 
( erfc(nz) / erfc(nR s ) for z > R s 



(25) 
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with the constant ( given by Eq. (|15"D. Clearly, in the limit R c — > oo (q — > 0), corresponding 
to a flat wall, Eq. (^) reduces to the previously derived result, Eq. (|14]). It is a remarkable 
feature that all effects of curvature are taken into account by the simple geometrical prefactor 
R c /(z + R c ), for sufficently small size ratios q. In this respect, the above result bears close 
similarity to the well-known Derjaguin approximation [25|. 



III. SIMULATION 

A. The simulation model 

In order to check the theoretical prediction of the forces at hard objects, we performed a 
monomer-resolved Molecular Dynamics (MD) simulation and calculated the mean force at 
the center of the star polymer to compare the data with theory. The model is based on the 



ideas of simulation methods applied on linear polymers and on a single star |31|. The 
main features are as follows. 

A purely repulsive and truncated Lennard- Jones potential acts between all Nf monomers 
at distances r: 



V LJ (r) 




m 12 -(^) 6 + l for r < 2^a LJ ; 

(26) 

for r > 2 1 /V LJ . 



Here, <tlj is the microscopic length scale of the beads and e sets the energy scale. In 
accordance with previous work ||16|| , we have chosen T = 1.2e//cs. 

An attractive FENE (finite extensible nonlinear elastic) potential additionally acts be- 



tween neighboring monomers along a chain |30 

I -15e(^' 2 

Q \ err t i 

^FENEV 



— 15e (^y) 2 hi 



Ro J 



for r < R ; 

(27) 

oo for r > R . 



This interaction diverges at r = R , which determines the maximal relative displacement of 
two neighboring beads. The energy e is the same as in Eq. (p6|), whereas for the length scale 
Ro we have chosen the value Ro = 1.5 olj. 

To accommodate the polymer arms, a hard core with radius Rd is introduced at the center 
of the star; its size depends on the arm number /. Accordingly, the interactions between 
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the monomers and the central particle were introduced. All monomers had a repulsive 
interaction V£j(r) of the truncated and shifted Lennard- Jones type with the central particle, 

{oo for r < R d ; 

(28) 
V LJ (r-R d ) forr>R i , 

whereas the innermost monomers in the chain an additional attractive potential V A p ENE (r) 
of the FENE type, namely 

oo for r < R d ; 

V F P mE (r) = { ~ (29) 

Vfene(^ - Rd) for r > R d . 

Finally, all monomers interact with the colloid or with the wall by a hard potential. We note 
that exactly this simulation model was already used by Grest et al. in their simulations of 
linear and star polymers in good solvent conditions [ROL |31|j. 



The timestep is typically At = 0.002t* with t* = y/ma^/e being the associated time unit 
and m the monomer mass. After a long equilibration time (500 000 MD steps), the mean 
force at the core of the star whose center is held at the position R and its dependence on 
the arm number / separations, is calculated as the expectation value over all instantaneous 
forces acting on the star core: 



F(R) = ( -V R 



' fN f 

E^jd^-RD+E^FENEd^-Rl 

k=l 1=1 



(30) 



where the first sum is carried over all fN monomers of the star and the second only over the / 
innermost monomers of its chains. The direct force between the central particle and the wall 
did not need to be considered, as the center-to-surface distance was always kept at values 
where this force was vanishingly small. Note that choosing the origin of the coordinate 
system on the surface of the colloidal particle or wall, at the point of nearest separation 
between the star center and this surface, and also the z-axis in the direction connecting this 
origin with the star center, we immediately obtain R = |R| = z. 

We have carried out simulations for a variety of arm numbers / and size ratios q, allowing 
us to make systematic predictions for the /- and g-dependencies of all theoretical parameters. 
In attempting to compare the simulation results with the theoretical predictions, one last 
obstacle must be removed: in theory, the fundamental length scale characterizing the star 



12 



is the corona radius R s . The latter, however, is not directly measurable in a simulation in 
which, instead, we can only assess to the radius of gyration R g . Yet, we have previously 
found that the ratio between the two remains fixed for all considered arm numbers /, having 
the value R s /R g ~ 0.66 [TE\. We now proceed with the presentation of our MD results. 



B. Star-wall and star-colloid interactions 



We consider at first a star polymer near a hard wall. The theoretical prediction of 



the effective interaction force is given in Eq. (0). First, we consider the limit of small 
separations, z — > 0, which allows us on the one hand to test the theoretical prediction 
F sw (z) = ksTAf 3 / 2 /z there and on the other hand to fix the value of the prefactor A, which is 
expected to have in general a weak /-dependence. For this prefactor, some semi-quantitative 
theoretical predictions already exist: For / = 1, 2 the prefactor may be calculated from the 
bulk and the ordinary surface critical exponents u, 7 and 7 , 7° of the n-vector model. For 
n = this results in A(/ = 1) = (7 - H)/v and 2 3 / 2 A(/ = 2) = (7 - 7 °)/v = l/v 



321, |33[. Numerical values for the exponents are known from renormalization group theory 



and simulation fl34| , p5|1 and yield A(/ = 1) ~ 0.83 and A(/ = 2) « 0.60. On the other hand, 
for very large functionalities, / 3> 1, one can make an analogy between a star at distance 
z from a wall and two star polymers whose centers are kept at distance r = 2z from each 



other p0[ . Indeed, for very large /, the conformations assumed by two stars brought close 
to each other is one in which the chains of each star retract to the half-space where the 
center of the star lies, a situation very similar to the star-wall case. Then, one can make the 
approximation F sw (z) = F ss (2z), where F ss denotes the star-star force. For the latter, it is 



known 121 that it has the form: 



F ss (r) = (r - 0), (31) 

implying for the coefficient A the asymptotic behavior: 

lim A(/) = Aoo = A = 0.14. (32) 
f^oo Sb 

Since there is no theory concerning the values of A in the intermediate regime of /, A is 
used as fit parameter. Its value can be obtained by plotting the inverse force 1 / F sw (z) against 
z for small separations z to the hard wall. The results are shown in Fig. |j. Looking first 
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at the inset, we see that, as for the earlier case of star-star interactions [fTE], the reciprocal 
force curves do not go through the origin, as a result of the finite core size, Rd- Once this is 
subtracted, though, straight lines passing through the origin are obtained, verifying in this 
way the 1 / z-behavior of the force and the associated logarithmic dependence of the effective 
potential at small separations. The values for A(/) can be immediately read off from the 
slope of the curves and they are summarized in Table |. There and in Fig. [5] we see that A is 
indeed a decreasing function of / but the asymptotic value Aoo = 5/36 is still not achieved 
at arm numbers as high as / = 100. 

The decay parameter k is fixed by looking at the force at larger separations and the 
obtained are also summarized in Table | and shown in Fig. [5|. As expected, k is of the order 
Rg 1 , as witnessed by the fact that the product kRs is of order unity. A monotonic increase 
of K&g with the arm number / is observed, consistent with the view that for large / stars 
form compact objects with an increasingly small diffuse layer beyond their coronae |12| . 



With parameters A and k once and for all fixed from the star-wall case, we now turn our 
attention to the interaction of a star polymer at a hard sphere of finite radius R c , equivalently 
size ratios q ^ 0. Here, the force is given by the full expressions of Eqs. (pi]), (E^) and (f23"D; 
for small enough size ratios q, the approximation fts max — > oo gives rise to a simplified 
expression for the force and to the analytical formula, Eq. fl25]) for the effective star-colloid 
potential. Our purpose is twofold: to test the validity of these simplified expressions as a 
function of q and also to find an economical way to parameterize s max as a function of q 
for those values of the size ratio for which the approximation Ks max — > oo turns out to be 
unsatisfactory. 

We show representative results for fixed arm number / = 18 and varying q in Fig. ||; 
results for different /-values are similar. It can be seen that the simplified result arising 
from allowing s max — > oo yields excellent results up to size ratios q < 0.3, see Figs. ||(a) 
and (b). However, above this value, the approximation of integrating the osmotic pressure 
up to infinitely large distances breaks down, as it produces effective forces that are larger 
than the simulation results, especially at distances z of order of the radius of gyration R g . 
These are the dashed lines shown in Figs, ^(c)-(e). The overestimation of the force is not 
surprising: as can be seen from Fig. ^ and Eq. (|19"D, we are integrating a positive quantity 
beyond the physically allowed limits and this will inadvertently enhance the resulting force. 
Hence, we have to impose a finite upper limit s max for size ratios q > 0.3 in order to truncate 
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the contribution of the Gaussian tail in the integral of the osmotic pressure in Eq. (|19D . 

In Fig. |7| a typical snapshot of a star polymer at a colloid illustrates the situation. One 
can see that the main contribution of the osmotic pressure results from in the inner region of 
the star. The outer region of the chains only interact weakly with the sphere. The question 
now is how the value of s max must be chosen. As can be seen from Eq. (|2D|), this quantity 
is dependent on, R s , z and S max . (The latter depending on q means that # max and q should 
not be treated as independent quantities.) It would be indeed most inconvenient if for every 
combination of these we would have to choose a different upper integration limit. Hence, 
we have attempted to transfer all dependence of s max onto the maximum integration angle 
#max- We found that this is indeed possible and, in fact, the angle max (q) has a very weak 
g-dependence: starting with a value # max ~ 45° at q = 0.3, we find that it then quickly 
saturates into the value # max ~ 30° for all q > 0.35. In this way, we are able to obtain the 
corrected curves denoted by the solid lines in Figs. ||(c)-(e), showing excellent agreement 
with the simulation results. 

We finally turn our attention to the /-dependence of the forces for a fixed value of the 
size ratio, q = 0.33. In Fig. [8] we show the simulation results compared with theory for 
a wide range of arm numbers, 5 < / < 50. For the theoretical fits, the values of A and 
k from Table | were used, whereas the value of the maximum integration angle was kept 
fixed at (9 max = 30° for all /-values. The agreement between theory and simulation is very 
satisfactory. 

Thus, our conclusions for the star polymer-colloid interaction read as follows: the general, 
analytical expression for the force between the two is given by Eqs. (|21~1) , ( p2|) and (p3|). 
supplemented by Eq. (p0|) in which the angle 9 max has to be chosen as discussed above for 
q > 0.3. An analytical formula for the effective interaction potential V sc (z) is not possible for 
such size ratios. Rather, the results for the effective force have to be integrated numerically 
in order to obtain V sc (z). For size ratios q < 0.3 on the other hand, the approximation 
•Smax — ► oo in Eqs. (plf), (|22| ) and (^3|) for the effective force can be made, thereby also 
allowing us to derive a simple, accurate, and analytic form for the interaction potential 
between a star polymer and a colloid, given by Eq. fl25|). These results form the basis of 
the statistical-mechanical treatment of star polymer-colloid mixtures in terms of standard 
liquid-state theories; the availability of analytical results for the pair interactions greatly 
facilitates the latter. A many-body theory of star polymer-colloid mixtures was put forward 
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recently by Dzubiella et al. [28[|, who employed the above-mentioned effective interactions 
in order to study the fluid-fluid separation (demixing transition) in such systems. The very 
good agreement with experimental results obtained in that work offers further corroboration 
of the validity of the interactions presented here. 

As the ultimate goal of the derivation of the interactions we present here is precisely 
to allow theoretical investigations of star polymer-colloid mixtures, we present in the next 
section a short account of a revision of the star-star interactions for the case of very low arm 
numbers. In this way, mixtures containing stars with arbitrary arm numbers, ranging from 
free chains (/ = 1, 2) to the "colloidal limit" of / ^> 1 can be studied in full generality. 



IV. REVISED STAR-STAR INTERACTION FOR SMALL ARM NUMBERS 



The effective interaction between two stars in a good solvent was recently derived by 
theoretical scaling arguments and verified by neutron scattering and molecular simulation 
|12| , [TJ|, [15], ITS] , leading thereafter to the phase diagram of the system |3(| |37J. The pair 
potential was modeled by an ultrasoft interaction which is logarithmic for an inner core and 



shows a Yukawa-type exponential decay at larger distances |12, |36 



for r < a s ; 



| N 1 ^A_ exp (_g( r _ as )) for r >a, 



(33) 



l+vT/2 

However, the theoretical approach giving rise to Eq. fl3"3]) does not hold for arm numbers 



/ < 10, because the Daoud-Cotton model of a star pT| , on which the Yukawa decay rests, is 
not valid for small /. In these cases, the interaction has to a shorter-ranged decay for r > a s . 
The shortcomings of the blob model can be made evident if one considers the extreme limit 
/ = 1, corresponding to free chains. There, the geometrical blob picture and the associated 



"cone approximation" |38[ break down. It is therefore instructive to consider known results 
about the effective interactions between free chains in order to obtain some insight for the 
case at hand. 

Most of the work done on chain-chain interactions concerns the effective potential between 



the centers of mass of the chains P, 10, 39, 40, 41]. Theoretical approaches considering two 



chains JO], simulations of two chains |39|, |40"fl, as well as recent, state-of-the-art simulations 



of many-chain systems j5L 10 1 all reach the conclusion that the effective center-of-mass to 
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center-of-mass interaction has a Gaussian form with its range set by the radius of gyration 
of the chains. Here, we are interested in a slightly different interaction, namely that between 
the end-monomer of one chain and the end-monomer of the other. However, at distances of 
the order of R g or larger, whether the centers of mass or the end-monomers choice of the two 
chains are held fixed should not make much difference. Therefore, we assume a Gaussian 
decay of the star-star potential for small /-values and center-to-center distances larger than 
a s . We emphasize that only the large distance decay of the interaction is affected; its form 
at close approaches has to remain logarithmic [^6|. Accordingly, we propose the following 
star-star pair potential for arm numbers / < 10, replacing the Yukawa by a Gaussian decay: 

5 , — ln(— ) + k-st for r < a s ; 



2^ exp (-r 2 (r 2 - a 2 )) for r > a s 



where r(/) is a free parameter of the order of 1/R g and is obtained by fitting to computer 
simulation results, see Fig. ^. For / = 2 we obtain the value r = 1.03 which, together with 
the potential in Eq. Q3~4] ) above yields for the second virial coefficient of polymer solutions the 
value B 2 /Rg = 5.59, in agreement with the estimate 5.5 < B 2 / R? g < 5.9 from renormalization 
group and simulations HlQ|] . For / = 5 we find r = 1.12, which leads to B 2 j R\ = 11.48, in 



accordance with Monte Carlo simulation results |42|, [43 . 

The very good agreement between the logarithmic-Gauss-potential of Eq. ( |34D and the 
simulation data for / < 10 can be seen in Fig. |9|. In the inset of this Figure, it can also be 
seen that the Yukawa decay is way too slow there. Hence, the potential of Eq. (|33|), has a 
longer range than the true interaction for small /, a property that explains the discrepancies 
between the simulated and theoretical second virial coefficients based on this potential, which 
have been reported by Rubio and Freire |EJ in their numerical study of low-functionality 
stars. At the same time, with increasing /, the roles of the Gaussian- and Yukawa-decays 
are reversed: in Fig. [H], we show simulation and theory results for / = 10. The original, 
logarithmic- Yukawa potential brings about better agreement now, as already established by 
earlier studies on stars with high arm numbers |T4], [Tffl. To summarize, we propose 
two analytic expressions for the effective star-star potential, valid in complementary regimes 
of the functionality /. The first one concerns the regime / < 10 with the validity of the 
logarithmic-Gauss-potential of Eq. (j3~4"D being established; in the second regime, / > 10, the 
logarithmic- Yukawa-potential of Eq. ( |3~3"D holds. We remark that the ultimative decay of the 
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effective interaction for very long distances is still Gaussian even for very large /, but this 
is not relevant for B2 as it occurs for much larger distances than the corona diameter. 

V. SUMMARY AND CONCLUDING REMARKS 

In summary, we have presented analytic results for the force between a colloid and a star 
polymer in a good solvent, accompanied with an analytic expression for the corresponding 
pair potential which is valid for size ratios q < 0.3. The validity of these expressions was 
established by direct comparison with Molecular Dynamics simulations. It should be noted 
that our theoretical approach is in principle generalizable to arbitrary geometrical shapes 
for the hard particle, thus opening up the possibility for studying effective forces between 
stars and hard ellipsoids, platelets etc. Further, a revised form for the star-star interaction 
for small functionalities has been presented, while at the same time the logarithmic- Yukawa 
form of this interaction remains valid for functionalities / > 10. 

The practical advantage of the present results is that they greatly facilitate the study 
of structural and thermodynamic properties of concentrated star polymer-colloid mixtures; 
a first attempt towards this has already been undertaken ||28|| . In this work, we limited 
ourselves to the case where the star is smaller than the colloid, i.e., q < 1. The study of 
the inverse case may be possible by applying the ideas presented here, however additional 
complications arise through the possibility of the star to "surround" the smaller, colloidal 
particle, in which case one part of the arms acts to bring about a repulsion with the colloid 
and another causes an effective attraction between the two. Furthermore, a pair poten- 
tial picture for the many-body system become more and more questionable for larger q as 
effective many-body forces |19j will play a more dominant role in this case. 
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/ 


Rd/R s 


A 


kR s 


2 


0.006 


0.46 


0.58 


5 


0.018 


0.35 


0.68 


10 


0.06 


0.30 


0.74 


15 


0.12 


0.28 


0.76 


18 


0.09 


0.27 


0.77 


30 


0.12 


0.24 


0.83 


40 


0.152 


0.24 


0.85 


50 


0.152 


0.23 


0.86 


80 


0.273 


0.22 


0.88 


100 


0.303 


0.22 


0.89 



TABLE I: The fit parameters arising from the comparison between theory and simula- 
tion for the star-wall and star-colloid interaction. The values of Ra shown here are not 
exactly the same as the input core size; they are just in the same order of magnitude, 
deviating only slightly from the real input value. They are still corresponding to mi- 
croscopic length, and are thus irrelevant at length scales r ~ a s . A is the overall pref- 
actor and k the inverse Gaussian decay length, both used in Eqs. (|2lj ) and (|25|) . a s = 
2R S = 0.66 cr g denotes the corona diameter of the stars, as measured during the simulation. 
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FIG. 1: Star polymer (black-shadowed particle) interacting with a flat wall. The star poly- 
mer consists of a inner core region, where the scaling behavior is dominant, whereas the 
outer regime is shadowed and indicates the exponential decay of the osmotic pressure. 
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FIG. 2: Snapshot of a simulation showing a star polymer interact- 
ing with a flat wall, at a small center-to-surface distance. The mirror- 
reflected image of the star, on the right, helps demonstrate that the con- 
figuration is similar to that of an isolated star with twice as many arms. 
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FIG. 3: Star polymer (black-shadowed particle) interacting with a colloidal particle (grey 
sphere). The dark and shadowed regions of the star have the same meaning as in Fig. |l|. 
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FIG. 4: Reciprocal effective force between a star polymer and a hard flat wall plotted against 
the distance z between the star center to the surface of the wall for small z-values. The 
dependence F(z) ~ \jz is confirmed by the simulation results (symbols). The prefactor 
of the potential depends on / and manifests itself in the different slopes of the reciprocal 
forces. The inserted plot shows the divergence of the force at the distance z = Rd, which is 
subtracted from z in the outset of the plot, to achieve divergence of the force in the origin. 
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FIG. 5: The pref actor of A and the decay parameter k of Eq. (f[3|) plotted against the 
functionality /. The value of A = 5/36 ~ 0.14 for / 3> 1 is not reached but the simulation 
data tend to this value very slowly, k shows a monotonic increase with arm number /. 
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FIG. 6: Comparison between simulation (symbols) and theoretical (lines) results for 
the effective force between a star polymer and a colloidal particle for different size ra- 
tios q, as a function of the center-to-surface separation z. The arm number here is 
/ = 18. The solid lines in (a) and (b) are derived from Eq. ([H]) for s max — > oo. 
In (c)-(e) the curves derived by means of this approximation are shown dashed and 
they increasingly deviate from the simulation results as q grows. Thereby, a finite up- 
per integration limit has to be introduced (see the text), producing the curves denoted 
by the solid lines in (c)-(e) and bringing about excellent agreement with simulation. 



29 



FIG. 7: Snapshot of a typical configuration of a star polymer with / = 18 arms near a 
colloidal sphere with q = 0.75. One should notice that predominantly the inner region of the 
star interacts with the hard sphere, yielding the main contribution of the inner core regime to 
the osmotic pressure of a region, determined by 6* max ~ 30°. Thereby, the upper integration 
limit s max in Eq. (|19|) is limited, see also the geometrical aspects of Eq. (|20|) and Fig. |3|. 
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FIG. 8: The effective force between a star polymer and a colloid for different arm numbers 
/ and q = 0.33 plotted against z, the distance of the star center to the surface of the colloid. 
The lines are the theoretical and the symbols the simulation results. For clarity, the data have 
been shifted upwards by constants: / = 10 : 1, / = 18 : 2, / = 30 : 3, / = 40 : 4, / = 50 : 5. 
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FIG. 9: Effective force between two star polymers plotted against the center-to-center dis- 
tance r for arm numbers / = 2 and / = 6. The simulation data (symbols) coincide with 
the logarithmic-Gauss expression (solid lines) of Eq. (j34[) . In the inset, the outer distance 
region is enlarged in order to clearly show the validity of the Gaussian decay in this /- 
regime (solid line), whereas the Yukawa form (dashed line) produces poor agreement there. 
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(r-2R d )/a g 

FIG. 10: Effective force between two star polymers plotted against the center-to- 
center distance r for arm number / = 10. In contrast to Fig. || here the 
Yukawa form (solid line) gives an accurate description of the decay of the in- 
teraction at large separations, whereas the Gaussian form (dashed line) does not. 
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